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Work under contract the NASA administered research contract NAS5-32073 to . 

Computer Sciences Corporation (CSC) was completed under a consultant arrangement f - 23 

between CSC and the Principle Investigator (PI), Dr. Chris R. Shrader (CSC employee from 
7/85 - 3/94). This report describes the final results of work performed by the PI on this task. 

Appended to the report are several resulting publications. The work described has been 
performed during weekend and evening hours on an IBM compatible 486/33 personal ! 

computer, utilizing the IDL Ver 3.5.1 for windows software purchased privately by the PI. r 
Also appended is an Invoice Detail in which the hours being charged to the task are tabulated*. 


The scientific focus of this program was the broad-band spectroscopic study of the 
subclass of quasi stellar objects known as "Blazars" using, primarily, data obtained with 
Compton Gamma-Ray Observatory (CGRO) and the International Ultraviolet Explorer (IUE), 
and ground-based optical and radio observatories. Of interest are episodic patterns of extreme 
variability, or outbursts in a number of sources which we monitor photometrically on a 
regular basis. This has led to coverage of a number of major events over the last several 
years. In several cases our observations and analysis have already been presented in scientific 
journals and or at scientific meetings. I include as an appendix to this document copies of 
several journal articles which have resulted. In one case involving the Blazar 3C 345, we 
were able to constrain Blazar emission models which incorporate an accretion flow onto a 
central compact object and a jet of relativistic particles. More fundamentally, several 
empirical trends in the multi-wavelength light curves were established - the UV luminosity 
associated with the accretion flow onto the central black hole was found to be correlated with 
the variability of the underlying synchrotron component, which is widely believed to be 
responsible for the radio to high-energy gamma-ray emission. Nominally, this suggests that 
Compton up-scattering of optical and UV photons may comprise a significant component of 
the high-energy emission. 


The most recent Blazar outburst, for which the data analysis is still ongoing involved 
the source known as OJ 287 (this event occurred during November 1994). We were able to 
obtain high-energy gamma-ray observations at the nominal peak of the outburst. This resulted 
in the detection of a source which was clearly undetectable at quiescence, thus providing new 
insight into the relationship between the high-energy emission and the lower-energy broad- 
band continuum - one of the "holy-grails" of current quasar research. What makes our 
observations even more intriguing, is that historically, this source has exhibited a very , 
approximate 10-year "periodicity" in its outburst behavior, based on which activity was 
expected for 1994. 


The data obtained for OJ 287 have thus far been analyzed in only a somewhat 
preliminary manner. It is expected however that the results will be fully analyzed and 
prepared for publication bv the summer 1995 time-frame. All other work on this project has 
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ABSTRACT 

An ongoing program of photometric monitoring of selected active galactic nuclei has revealed a 
significant outburst in the optically violently variable QSO 3C 279. The event occurred during May and 
June of 1992, and was characterized by a peak amplitude of 2.3 mag in the R band. We obtained nearly 
simultaneous optical and ultraviolet spectroscopy from 1 150 to 8600 A of the source at the approximate 
midpoint of the outburst. The outburst spectrum is well represented by a power law with index 
a. = — 1.77 (where F v cc v a ) which is relatively flat in comparison with typical quiescent spectra for this 
source. No emission features are discemable in our data, nor is there any significant Lyman-edge 
absorption. The observations and data analysis procedures are described, and the resulting photometric 
light curves and spectra are presented. We also present a long baseline photometric light curve for 
context. Comparisons with previously published optical-UV outburst and quiescent spectra are made, 
and luminosity-dependent spectral variations are discussed. 


1. INTRODUCTION 

The QSO 3C 279 is a well known member of the class of 
optically-violent-variable quasars, or OWs which are a 
subclass of the so-called blazars. It is among the most dra- 
matic examples of the blazar class having a well docu- 
mented history of optical outbursts, including flares of ~ 5 
mag lasting up to several years and involving energy out- 
puts of — 10 55 ergs. (Webb et al. 1990). 3C 279 is perhaps 
the prototypical superluminal radio source (MofiFet et al. 
1972; Unwin et al. 1989). Additionally it is now well 
known for being the first quasar detected at energies 
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greater than 1 GeV by the EGRET experiment on the 
Compton Gamma Ray Observatory. During one observa- 
tion it exhibited the highest apparent brightness of any 
quasar yet detected by EGRET and was found to be highly 
variable at high energies on time scales as short as a few 
days and as long as (several) months (Hartman et al. 
1991; Kniffen et al. 1993). Unfortunately there is as yet no 
well established general correlation between flaring in the 
IR-UV and the global continuum properties for this or any 
other blazar, nor are the broadband characteristics of low 
amplitude variability particularly well established. 
Progress has been made in attempts to extend current the- 
oretical models, which have been reasonably successful in 
accounting for the radio-UV spectra to the ~l-20 keV 
x-ray domain, and more recently into the —GeV domain 
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3C 279 Photometric Light Curve 

- — : — — ’ — i — ■ — ■ — ■ — i — ■ — . — . — r— — ■ — 1 — i — ■- 


1 1 

: — ,■ * ■ — v — ; 1 »• 1 1 — 1 1 1 : 


> t C_ «- 


: ^ ° “ 

12 

- S 0. f -. 

1 1 1 1 


; A *= R-bond ; 

13 

_ X = V-bond j 

A 3 

CJ 

■o 

4 4 

3 

L ** J 

CP 

o 

> XX * * ] 

2 

AX : 

15 

* - 


A A 


A : 

16 

F , ,, 4 , j 


iSS 4 : 

17 

L_. — . — . — i — . — . — « • • ■ i . , 


8720 8740 8760 8780 8800 8820 8840 

JD - 2440000 


Fig. 1. R and V band photometric light curves for Spring 1992. Epochs 
for which optical and UV spectroscopic measurements were obtained are 
indicated. 

(e.g., Marscher et al. 1992; Dermer et al. 1992). However, 
all such efforts are hampered by the current lack of a clear 
empirical multiwavelength variability picture. The dra- 
matic outbursts seen in the IR-optical-UV are likely to be 
a central element in any comprehensive understanding of 
the blazar phenomenon, given the remarkable energetics 
involved coupled with the rapid, and thus observable, ev- 
olution of these events. 

We have established a program of optical photometric 
monitoring, which is described in Sec. 2. This program led 
to the discovery of an outburst of 3C 279 during the spring 
of 1992. Our observations of the source revealed a bright- 
ening by several magnitudes from late March to early June 
1992 (JD 2448740-8817; Fig. 1), which in turn triggered 
our target-of-opportunity spectroscopic observations in the 
optical and UV. The March, preoutburst level represents 
the nominal quiescent level for this source. The brightening 
of = 2.3 mag can be characterized as an event of interme- 
diate scale, referring to the RHO long baseline light curve 
in Fig. 2 for an historical perspective. The mean AR/At 
over the ten day interval, JD 2448760 : 8770, which shortly 
preceded our spectroscopic observations was —0.13 mag/ 
day. 

In Sec. 2 we describe our photometric and spectroscopic 
observing programs and detail our coverage of the 1992 
outburst of 3C 279. In Sec. 3 we describe details of our 
methods of data analysis. Additionally, we offer interpre- 
tation of our data and discuss the general optical-UV spec- 
tral dependence on overall luminosity. 


3C 279 Long Term Light Curve 



JD - 2400000 


Fig. 2. Long term B band photometric light curve. The 1992 flare is 
labeled (lower left). 

2. OBSERVATIONS 
2.1 Photometry 

Our program of photometric monitoring utilizes several 
small ground-based telescopes located at the Rosemary 
Hill Observatory operated by the University of Florida 
(RHO) and the Foggy Bottom Observatory (FBO) at 
Colgate University. Monitoring by these instruments al- 
lows us to promptly identify periods of activity in an object 
among our blazar sample, which has been preselected on 
the basis of average brightness and variability history. This 
information then allows us to trigger target-of-opportunity 
programs of followup spectroscopy in the optical and UV 
(Table 1). 

We obtained optical photometry at FBO using the 16 
inch telescope and Photometries' PM3000 CCD system. 
Our measurements are summarized in Table 2. Each value 
represents the average of, typically, three to five exposures 
of three minutes duration. No magnitude-transformation 
corrections have been made since they are negligibly small. 
For this telescope and detector system, the transformation 
coefficients are as follows: 

R — r= 0.00 ±0.02 ( K—R), 

V— R = l. 04 ±0.02(v—r). 

The comparison star used, located about 02 '52" south and 
< 14" west of the quasar ( V— 12.42), has a color index of 
V—R =0.365, where V and R refer, respectively, to the 
Johnson and Cousins photometric systems. Since the qua- 


Table 1. Log of spectroscopic observations. 


Observatory 

Instrument 

Bandpass/resol ution 

Exposure time 

GMT (1992)/JD 
epoch 

IUE 

LWP 

AA1950-3350/7 A 

180 min 

May 29.80/8772.3 


SWP 

AA 11 50-2000/6 A 

180 min 

May 29.67/8772.2 

McDonald 

Cassegrain 

spectrograph 

AA7138-8600/5 A 
AA4620-7380/10 A 

10 min 
20 min 

June 1.25/8774.8 
June 3.07/8776.7 
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. Table 2. Columns are Julian Day —244000, magnitude, uncertainty, 
photometric band, and observatory ID. 


JD 

mag 

o 

fil 

Obs 

8665.7 

16.185 

0.152 

R 

FBO 

8685.8 

14.665 

0.028 

R 

FBO 

8703.6 

15.171 

0.052 

R 

FBO 

8704.5 

15.120 

0.170 

V 

RHO 

8706.7 

14.834 

0.026 

R 

FBO 

8711.6 

15.080 

0.028 

R 

FBO 

8716.7 

15.346 

0.031 

R 

FBO 

8717.7 

15.408 

0.031 

R 

FBO 

8718.7 

15.502 

0.034 

R 

FBO 

8737.5 

16.400 

0.060 

V 

RHO 

8740.6 

15.484 

0.032 

R 

FBO 

8741.7 

15.309 

0.030 

R 

FBO 

8749.6 

15.295 

0.039 

R 

FBO 

8753.5 

15.120 

0.041 

R 

FBO 

8761.5 

14.742 

0.023 

R 

FBO 

8762.5 

15.300 

0.080 

V 

RHO 

8762.6 

14.477 

0.022 

R 

FBO 

8763.5 

14.233 

0.021 

R 

FBO 

8763.6 

14.730 

0.013 

V 

RHO 

8764.6 

13.868 

0.020 

R 

FBO 

8764.6 

14.388 

0.012 

V 

RHO 

8765.5 

13.886 

0.021 

R 

FBO 

8765.6 

14.357 

0.012 

V 

RHO 

8767.5 

14.520 

0.100 

V 

RHO 

8768.5 

13.738 

0.021 

R 

FBO 

8768.5 

14.590 

0.080 

V 

RHO 

8768.6 

14.280 

0.019 

V 

RHO 

8771.6 

13.760 

0.022 

R 

FBO 

8771.6 

14.279 

0.019 

V 

RHO 

8772.5 

13.785 

0.021 

R 

FBO 

8776.5 

-13.782 

0.020 

R 

FBO 

8776.5 

14.344 

0.027 

V 

RHO 

8777.5 

13.915 

0.021 

R 

FBO 

8777.5 

14.418 

0.012 

V 

RHO 

8785.5 

13.413 

0.022 

R 

FBO 

8785.5 

13.864 

0.056 

V 

RHO 

8786.5 

13.064 

0.023 

R 

FBO 

8786.5 

13.529 

0.052 

V 

RHO 

8788.5 

13.705 

0.047 

R 

FBO 

8789.5 

13.369 

0.030 

R 

FBO 

8789.5 

13.899 

0.016 

y 

RHO 

8790.5 

13.653 

0.021 

R 

FBO 

8791.5 

14.330 

0.054 

R 

FBO 

8791.5 

14.930 

0.100 

V 

RHO 

8806.5 

15.950 

0.102 

R 

FBO 

8815.5 

16.200 

0.201 

R 

FBO 


sar has a v—r color index between 0.45 and 0.55 during 
the period of our observations, the color difference was 
A( F— i?)~0.13±0.05. The color correction for V—v 
would thus be approximately 0.005 mag, which is negligi- 
ble in comparison with our typical la measurement error. 
The quoted uncertainties for the FBO measurements in 
Table 2 include internal-statistical error resulting from the 
differential photometry and an uncertainty of 0.02 mag in 
the magnitude of the comparison star. 

Additional photometric coverage was obtained at RHO, 
where 3C 279 has been routinely monitored since May of 
1971 (Webb 1990; Pica et al. 1988). The RHO observa- 
tions are made in the international photographic system 
(PG) and are converted to Johnson B magnitudes using 


an empirically determined conversion factor of 0.3 mag 
(£— PG=0.3). 

2.2 Ultraviolet Spectroscopy 

We obtained UV spectra during the 1992 outburst at the 
epochs labeled on Fig. 1 . The UV observations were made 
with the IUE short and long wavelength instruments in 
low dispersion mode [see e.g., Newmark et al. (1992) for 
details on the IUE spacecraft and instrumentation]. These 
observations were implemented in a target-of-opportunity 
mode, intended specifically to obtain time-critical outburst 
observations of blazar AGN. 

We note that the source was too faint to be detected 
with the IUE Fine Error Sensor, thus no precisely simul- 
taneous photometric measurement was possible and the 
target was acquired using a blind-offset procedure. Com- 
parison between the empirically determined guide star po- 
sition and the calculated position suggested that the offset 
was accurate to within several arcseconds. The resulting 
spectra were both found to be well centered with respect to 
the major axis of the 10" X20" oval aperture, which is the 
direction perpendicular to the dispersion, thus we believe 
that the absolute flux calibration should be valid to within 
nominal limits (~10%). 

Initial processing of the data was performed by the 
GSFC staff using the standard IUE SIPS software (e.g., 
Newmark et al. 1992). We then used the software at the 
GSFC Regional Data Analysis Facility to computationally 
apply a Gaussian weighted extraction slit to the two- 
dimensional “line-by-line” images, resulting in flux and 
wavelength calibrated spectra. In addition to this Gaussian 
extraction process, spectra obtained by the standard IUE 
SIPS extraction procedure were analyzed for comparison. 
No gross differences in the continuum flux levels were dis- 
cemable. The SWP was somewhat under exposed, while 
the LWP exposure level was nominal. The observations are 
summarized in Table 1. 

2.3 Optical Spectroscopy 

Optical spectra were subsequently obtained using the 
Large Cassegrain Spectrograph and the CRAF Cassini 
1024 X 1024 CCD on the 2.7 m telescope of McDonald 
Observatory. On UT 1992 June 1 and 3 we covered the 
wavelength bands AA7 13 8-8600 and AA4620-7380 at spec- 
tral resolutions of 5 and 10 A, respectively. The red spec- 
trum is of lower signal to noise than the blue spectrum, as 
a result of the shorter exposure time and higher dispersion. 

Table 1 is a log of our spectroscopic observations. We 
note that the UV and optical spectroscopic observations 
were not strictly simultaneous, but were obtained within a 
4 day interval. Photometric coverage (Fig. 1), however, 
indicates that the source remained nearly constant during 
the interval. 

3. DATA ANALYSIS AND DISCUSSION 

Referring to the long baseline light curve (Fig. 2) and 
the spring 1992 light curve (Fig. 1) defining the event, it is 
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seen that the source was near its “quiescent” magnitude as 
of mid-March 1992, followed by a gradual increase in 
brightness of AR-0.5 mag over a two week span. It then 
brightened sharply, approximately 1.3 mag over about 
seven days. This increase was followed by a flattening of 
the light curve which lasted for some three weeks. It was 
during this “plateau” that we obtained optical and UV 
spectra. The mean AR/At over a ten day interval shortly 
preceding our spectroscopic observations was —0.13 mag/ 
day, which is small compared to the total amplitude of the 
event. A dramatic drop in brightness occurs subsequent to 
JD 2 448 790. Unfortunately, our sampling is sparse the 
during the majority of this decline, due to a combination of 
bad weather and impinging moonlight. Referring to Figs. 1 
and 2, it is interesting to note that the historical JD 242850 
(1938) outburst was followed approximately three years 
later by a moderate intensity peak comparable to the event 
we are reporting on here. Although this is most likely ran- 
dom coincidence, it is interesting to speculate that the vari- 
ability patterns are due to some stable (on a time scale of 
decades) underlying physical process, in which case we 
may see additional outburst activity in the 1995 time 
frame. 

A discontinuity of about 7% was seen in the region of 
overlap of our blue (AA4620-7380) and red (AA7138- 
8600) spectra, the latter being lower than the former. We 
referred the flux density scale to only one standard star on 
each night, and the nights appeared to be photometric. 
Generally, we expect the uncertainty in absolute flux den- 
sity to be a few percent, certainly less than 5%, so we 
interpret the disagreement in the overlapping region as 
real variability. For subsequent analysis, we have cor- 
rected the data for galactic extinction, using a color excess 
ofE(B-V) =0.02, and applying the correction algorithm 
of Seaton (1979). Additionally, we have applied a scaling 
factor to eliminate the 7% discontinuity and force the two 
optical spectra to connect smoothly. Regions containing 
atmospheric absorption bands at 7160, 7620, and 8200 A 
were excluded in subsequent analysis as well. We detect no 
statistically significant emission lines, noting that our data 
are signal to noise limited in the UV and red bands. At a 
redshift of z=0.538, Lya is well situated on a sensitive 
portion of the SWP detector, but the spectrum is poorly 
exposed and thus of low signal to noise. Based on our 
estimate of the signal to noise per resolution element, we 
estimate a la upper limit of o ew ~13 A. Civ (2.1550), 
typically the second strongest UV line seen in AGN spec- 
tra, is redshifted to a region of poor sensitivity on the 
LWP, however the level of the LWP exposure was more 
satisfactory than for the SWP. We estimate an upper limit 
of tr ew ~9 A for C iv emission. Additionally, we derived an 
upper limit of a ew =;2.1 A for H/3. We also examined each 
of the IUE archival spectra for 3C 279 and found no emis- 
sion lines. It is likely that the line-to-continuum ratio de- 
creases with increasing overall source brightness, however, 
the source in quiescence is generally undetectable by IUE. 
This lineless continuum is consistent with a synchrotron 
origin for the outburst, probably dominated by relativisti- 
cally beamed emission. We note, however that Webb et al. 


3C 279 Opticol/UV Light Curves 



Fig. 3. Recent photometric ( V band) and UV continuum (1800 A) 
history of 3C 279. The UV fluxes are in units of 10-14 erg/cm 2 /s/A. 

(1993) report no statistically significant line-continuum 
dependence on overall brightness for another OVV-QSO 
3C 345. 

Another feature characteristic of QSO spectra is an ex- 
cess above a power law shortward of ~5000 A, the so 
called “blue bump” (e.g., Malkan 1983), which for mod- 
erately red-shifted objects such as this, should be readily 
apparent in the blue to UV. Although this feature is typi- 
cally less pronounced in OVV quasars than in non-OVV 
quasars, it is certainly seen in some cases (e.g., Webb et al. 
1993; Malkan & Moore 1986). If present, this feature 
should be characterized by a flattening in the v-/v spec- 
trum plot in the vicinity of log(v) cr 14.8. No such feature 
is apparent in our data. To test this in a more quantitative 
manner, we fit two separate power laws to our composite 
spectra in v-/v space; one blueward and one redward of 
log(v) = 14.8. We found that at slightly better than the 4 a 
level the blue power law was steeper than the red power 
law. We conclude that no blue-bump component is present, 
the optical-UV light output being dominated by the blazar 
component. 


3C 279 Opticol-UV Continuum 5/29-6/1/92 



Fig. 4. Power law fit to the optical-UV continuum. The data are binned, 
transformed to the source frame, and corrected as described in Sec. 3. 
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Table 3. Uncorrected (col. 2-3) and corrected (col. 4-5) continuum 
fluxes. 


V 

F v 

a 

Ff' 


15.363 

-26.405 

0.085 

-26.382 

0.114 

15.301 

-26.266 

0.101 

-26.178 

0.114 

15.273 

-26.124 

0.086 

-26.099 

0.104 

15.247 

-26.106 

0.091 

-26.043 

0.101 

15.198 

-25.956 

0.064 

-25.953 

0.099 

15.135 

-25.971 

0.104 

-25.870 

0.105 

15.079 

-25.845 

0.064 

-25.795 

0.078 

15.046 

-25.763 

0.061 

-25.730 

0.066 

15.015 

-25.690 

0.063 

-25.642 

0.063 

14.979 

-25.579 

0.074 

-25.573 

0.074 

14.811 

-25.308 

0.011 

-25.272 

0.030 

14.805 

-25.286 

0.032 

-25.261 

0.032 

14.796 

-25.263 

0.034 

-25.243 

0.038 

14.787 

-25.253 

0.024 

-25.227 

0.027 

14.778 

-25.242 

0.021 

-25.211 

0.025 

14.770 

-25.226 

0.025 

-25.197 

0.025 

14.761 

-25.208 

0.022 

-25.183 

0.023 

14.753 

-25.197 

0.018 

-25.170 

0.018 

14.745 

-25.182 

0.015 

-25.158 

0.016 

14.737 

-25.176 

0.014 

-25.145 

0.018 

14.729 

-25.158 

0.014 

-25.133 

0.015 

14.721 

-25.143 

0.012 

-25.121 

0.013 

14.714 

-25.130 

0.013 

-25.109 

0.015 

14.706 

-25.116 

0.012 

-25.098 

0.016 

14.699 

-25.108 

0.011 

-25.086 

0.013 

14.692 

-25.098 

0.013 

-25.074 

0.014 

14.685 

-25.084 

0.012 

-25.063 

0.016 

14.678 

-25.079 

0.014 

-25.052 

0.017 

14.671 

-25.060 

0.011 

-25.041 

0.011 

14.664 

-25.051 

0.012 

-25.029 

0.012 

14.658 

-25.035 

0.012 

-25.018 

0.014 

14.651 

-25.030 

0.015 

-25.007 

0.016 

14.645 

-25.014 

0.012 

-24.995 

0.016 

14.638 

-24.998 

0.019 

-24.983 

0.019 

14.632 

-24.990 

0.018 

-24.971 

0.018 

14.626 

-24.977 

0.017 

-24.957 

0.018 

14.620 

-24.955 

0.022 

-24.943 

0.029 

14.614 

-25.003 

0.063 

-24.943 

0.063 

14.614 

-24.953 

0.025 

-24.943 

0.051 

14.608 

-24.982 

0.065 

-24.938 

0.065 

14.602 

-24.959 

0.061 

-24.921 

0.067 

14.596 

-24.969 

0.079 

-24.903 

0.079 

14.591 

-24.946 

0.066 

-24.889 

0.067 

14.585 

-24.915 

0.068 

-24.880 

0.069 

14.579 

-24.914 

0.065 

-24.874 

0.068 

14.574 

-24.913 

0.068 

-24.872 

0.069 

14.569 

-24.921 

0.075 

-24.872 

0.076 

14.563 

-24.937 

0.081 

-24.874 

0.082 

14.558 

-24.920 

0.079 

-24.874 

0.081 

14.553 

-24.899 

0.082 

-24.871 

0.084 


W also searched for the presence of a Lyman disconti- 
nuity in the UV spectrum, which would occur at 1402 A, 
but no statistically significant discontinuity was seen at 
that position in the spectrum. This implies that the cover- 
ing fraction of Lyman-edge absorbing material is zero to 
within measurement errors. 

Based on available data, there appears to be a positive 
correlation between the optical and UV flux in 3C 279; 
however, the UV measurement we obtained during the 
1992 outburst seems to be contrary to this trend. Figure 3 
shows the recent UV continuum (A obs ~ 1800 A) and pho- 


Power Low Index vs F 4400 



Fig. 5. Fv power law indices vs 4400 A flux per 10-14 ergs/cm 2 /s/A; 
X=IUE+ photometry, A = UBVRI, and (> = these observations. 


tometric ( F-band) history. Our data exhibit the lowest UV 
flux level among the data shown, despite the clear, sharp 
upturn in the photometric light curve. The JUE observa- 
tions took place approximately one week subsequent to the 
initial sharp rise in the optical light curve, so the possibility 
of a rise-time delay between the red and ultraviolet contin- 
uua cannot be ruled out. 

Our composite optical-UV spectrum is fitted by a power 
law of a = — 1.77 ±0.03 (where /va v°; Fig. 4). The spec- 
tral data are tabulated in Table 3, where we have included 
flux values corrected for galactic extinction, and for the 
cross-calibration scaling factor between the blue and red 
spectra as described above, as well as uncorrected fluxes. 
All of the data in Table 3 refer to the observers frame of 
reference. The power law index of 1.77 is significantly flat- 
ter than that seen in several cases where the source was 
observed in quiescence, but similar to the 1988 outburst 
cases reported for example by Makino et al. (1990) and 
Peterson et al. (1988). We calculated the source luminos- 
ity over the AA 1200-8600 A bandpass to be 5.4 X 10 46 
ergs/s for H 0 — 55 and g 0 = 0. This was obtained using the 
corrected flux values and assuming a smooth interpolation 
of the power law continuum distribution between the nom- 
inal 3250 A cutoff of the IUE LWP and the 4620 A blue 
cutoff of our optical spectrum. 

In an attempt to establish a relationship between bright- 
ness and optical-UV spectral index we have analyzed the 
available UV archival data, and identified cases where con- 
temporaneous photometric measurement are available 
from our monitoring database or from the published liter- 
ature. Additionally, we derived power law indices from 
published UBVRI photometric measurements (Webb et al. 
1990; Mead et al. 1990; Kidger et al. 1992). Figure 5 rep- 
resents a compilation of power law indices a (where 
fv ccv a ) versus the 4400 A continuum fluxes derived in 
this manner. If any trend at all is evident in Fig. 5, it is a 
flattening with increasing flux. Landau et al. (1986) pre- 
sented a similar plot where again there was no compelling 
evidence spectral-luminosity dependence. Bonnell et al. 
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(1993) suggest a flattening with increasing luminosity 
based on the UV data alone. 

We wish to acknowledge the support provided by ob- 
servatory staffs at IUE, RHO, and FBO. Parts of this work 


were supported under grants from the NASA IUE and 
Compton Observatory Guest Investigator Programs. S.G. 
participated as a visiting summer student at Colgate Uni- 
versity supported by a grant from the W. M. Keck Foun- 
dation to the Keck Northeast Astronomy Consortium. 


REFERENCES 


Bonnell, J. T., Vestrand, W. T., & Stacy, J. G. 1993, ApJ, 420 (in press) 
Dermer, C. D., Schlickeiser, R. & Mastichiadis, A. 1992, A.&A, 257, L27 
Hartman, R., el al. 1991, ApJ, 385, LI 

Kidger, M., Garcia Lario, P., & de Diego, J. A. 1992, A&AS 93, 39! 
Kniffen, el al. 1993, ApJ, 411, 133 
Landau, R., el al. 1986, ApJ, 308, 78 
Makino, el al. 1989, ApJ, 347, L9 
Malkan, M. A. 1983, ApJ, 268, 582 
Malkan, M. A., & Moore, R. L. 1986, ApJ, 300, 216 
Mead, A. R. G., Ballard, K. R„ Brand, P. W. J. L., Hough, J. H., Brindle, 
C., & Bailey, J. A. 1990, A&AS, 83, 183 
Moffet, A. T., Gubbay, J., Robertson, D. S., & Legg, A..J. 1972, in 
External Galaxies and Quasi-Stellar Objects, IAU Symposium No. 44, 
edited by D. S. Evans (Reidel, Dordrecht), p. 28 


Newmark, J. S., Holm, A. V., Imhoff, C. L., Oliversen, N. A., Pitts, R. E., 
& Sonnebom, G. 1992, NASA IUE Newslett. 47, p. 4 
Peterson, B. M., Wagner, R. M., & Korista, K. T. 1988 IAU Circ. No. 
4643 

Pica, A. J., Smith, A. G., Webb, J. R., Leacock, R. J., Clements, S., & 
Gombola, P. P. 1988, AJ, 96, 1215 
Seaton, M. 3. 1979, MNRAS 187, 73P 
Sitko, M., & Sitko, A. 1991, PASP, 103, 160 

Unwin, S. C., Cohen, M. H., Biretta, J. A., Hodges, M. W., & Zensus, J. 
A., 1989, ApJ, 340, 117 

Webb, J. R., Carini, M. T., Clements, S., Fajardo, S., Gombola, P. P., 
Leacock, R. J., Sadun, A. C., & Smith, A. G. 1990, AJ, 100, 1452 
Webb, J. R., ei al. 1993, ApJ (in press) 



The Astrophysical Journal, 422:570-585, 1994 February 20 

© 1994. Tbe American Aslronomical Society. All rights reserved. Printed in U.S.A. 


THE MULTIFREQUENCY SPECTRAL EVOLUTION OF BLAZAR 3C 345 
DURING THE 1991 OUTBURST 

J. R. Webb, 1 C. R. Shrader, 2 T. J. Balonek, 3 D. M. Crenshaw, 4 D. Kazanas, 2 S. Clements, 5 A. G. Smith, 5 
A. D. Nair, 5 R, J. Leacock, 5 P. P. Gombola, 5 A. Sadun, 6 H. R. Miller, 1 I. Robson, 8 R. Fujimoto, 9 
F. Makino, 9 T. Kii, 9 H, Aller, 10 M. Aller, 10 P. Hughes, 10 E. Valtaoja, 11 H. Terasranta, 1 1 
E. Salonen, 11 M. Tornikoski, 11 and W. Chism 1 
Received 1993 July 28; accepted 1993 August 30 

ABSTRACT 

The blazar 3C 345 underwent a 2.5 mag optical outburst between 1990 November and 1991 May. We have 
obtained 10 nearly simultaneous multifrequency spectra during the course of the outburst in order to study 
the multifrequency spectral variations of 3C 345 as a function of time. Although our observations were not 
sampled frequently enough to completely resolve the variations in every frequency band, the general rise and 
decline of the outburst were seen in the UV through radio with differing rise times. Simulations of an electron 
distribution injected into a tangled magnetic field show a relationship between frequency and characteristic 
timescale that was also observed in the radio variations of 3C 345. The two X-ray observations made during 
the monitoring period showed no evidence of variability. The multifrequency spectrum was modeled with two 
major components: a relativistic jet and a relativistic thermal accretion disk. Models calculated for each spec- 
trum indicate that the outburst can be explained in terms of these models by varying the high-energy cutoff of 
the injected electron distribution in the jet model, while also varying the mass accretion rate in the disk 
model. There is marginal evidence that the inferred accretion rate varies with the jet luminosity. 

Subject headings: quasars: individual (3C 345) — radio continuum: galaxies — - ultraviolet: galaxies — 

X-rays: galaxies 


1. INTRODUCTION 

The blazar 3C 345 ( z - 0.595; Burbidge 1965) has a long 
history of optical variations (Smyth & Wolstencroft 1970; 
Webb et al. 1988) and exhibits a high degree of variable polar- 
ization (Smith et al. 1986). Radio flux monitoring has shown 
that large-amplitude variations are common in that frequency 
range as well (Aller et al. 1985). In addition, superluminal radio 
components have been observed in 3C 345 (Kollgaard, Wardle, 
& Roberts 1989) which appear to be correlated with the radio 
flux variations. The general multifrequency spectrum of this 
source was studied by Bregman et al. (1986) and the variability 
of the radio to UV continuum was investigated by Brown et al. 
(1989). The previous studies found that the spectrum exhibits a 
“blue bump” in the optical-UV region and that the optical 
and IR spectral slope varies between — 1.5 and — 1.7. 

Prior to our observations, 3C 345 had been in an extremely 
faint state with a B magnitude of approximately 18.1 in 1990 
May (Kidger & Takalo 1990). In 1990 November, 3C 345 
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underwent an outburst brightening to nearly B = 15.40. Figure 
1 shows the long-term light curve of 3C 345 with the 1991 
outburst included. Except for some of the 1991 points, all of the 
observations in Figure 1 are photographic photometry from 
the University of Florida’s Rosemary Hill Observatory (RHO). 
Although the present outburst was not the brightest ever 
recorded for 3C 345, its relatively rapid increase from a 
minimum state and its subsequent rapid decline are remi- 
niscent of outbursts observed in this source during the early 
1970s. The 1991 outburst was followed by a second outburst 
which occurred in early 1992 (Kidger et al. 1993). 

The observations reported here were in a large part driven 
by an International Ultraviolet Explorer ( IUE ) Target of 
Opportunity (TOO) program designed to observe a quasar in 
outburst and to follow its evolution in the optical and UV. The 
TOO program focused on a small list of 15 optically violently 
variable (OVV) quasars and BL Lac objects which have shown 
large-amplitude optical flares during the years of monitoring at 
RHO. The monitoring frequency of the sources on this list 
were increased at RHO and at Colgate University's Foggy 
Bottom Observatory (FBO) in order to detect outbursts 
quickly. We activated the IUE TOO program based on obser- 
vations of3C 345 reported by Borgeest & Schramm (1991) and 
on follow-up photometry from RHO and FBO. The optical 
monitoring at RHO and FBO of 3C 345 was further increased 
to utilize every available clear night during the outburst. The 
X-ray observations were obtained by three of us (R. F., F. M., 
and T. K.) as part of a program of Ginga active galactic nucleus 
(AGN) observations. The radio, millimeter, IR, and supple- 
mentary optical observations were obtained under separate 
AGN monitoring programs and combined with our optical 
and UV data to form a comprehensive picture of 3C 345 
during this period of activity. 
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Fig. 1. — Long-term B-band curve of 3C 345 from Rosemary Hill Observatory. The plus signs along the top indicate the first day of each year. The duration of the 
1991 outburst is indicated by the bracket below the light curve. 


2. OBSERVATIONS 

In this section, we give a brief description of the observations 
in decreasing frequency order. References to data collection 
and analysis procedures in each frequency band are made 
where appropriate. A log of all of the observations is presented 
in Table 1. Column (1) lists the modified Julian Date, column 
(2) the civil date, column (3) the observatory, and column (4) 
the bandpass. 

2.1. X-Ray 

The X-ray observations of 3C 345 were obtained on 1991 
April 6.3— 6.5 UT, May 25.2-25.8 UT, and May 27.2-27.8 UT, 
with the large-area counter (LAC) array on Ginga (Turner et al. 
1989). The April 6.3-6.5 (UT) observation was fitted to a 
simple model consisting of a power law and a low-energy expo- 
nential absorption component, which yielded a photon index 
of 1.54 + 0.02 and a normalization of 2.03 x 10" 3 photons s' 1 
cm' 2 keV _1 at 1 keV. The May 25.2-25.8 and May 27.2-27.6 
observations yielded photon indices of 1.58 ± 0.05 and 
1.46 ± 0.23, and normalizations of 1.86 x 10~ 3 and 
1.27 x 10 3 between 1.16 keV and 23.16 keV. The individual 
May observations were combined; the weighted average of the 
photon index was 1.53 ± 0.05, and the normalization was 
1.589 x 10~ 3 . The hydrogen column density was fixed at 
1.0 x 10 20 cm' 2 for each observation, since the LAC in insen- 
sitive to column densities below 5 x 10 20 cm' 2 and since pre- 
vious estimations of the column density toward 3C 345 made 
with Einstein imaging proportional counter (IPC) observations 
(Webb 1988) also agree with Galactic column densities of 
around 1 x 10 20 cm' 2 . These values for N H are each consis- 
tent to within a nominal factor of 2 of the value obtained by 


applying an empirical scaling factor to the color excess 
E(B— V) = 0.01. 

2.2. Ultraviolet 

Ultraviolet observations were obtained with the Internation- 
al Ultraviolet Explorer (IUE) throughout the decline of the 
outburst. During a 200 day period, we obtained 12 good- 
quality low-dispersion spectra (six SWP and six LWP). The 
spectra were systematically reduced using the IUE RDAF 
Gaussian extraction routine (GEX). The most prominent fea- 
tures visible in the UV spectrum are the redshifted Lyman-a 
(1946 A) and the C iv line at 2470 A. The poor signal-to-noise 
ratio, especially in the SWP, and the presence of camera arti- 
facts and reseaux prevent firm identification of other lines. The 
identification of Lyman-/? 4- O hi in the 1620 A region by 
Bregman et al. (1986) is questionable, since there is a strong 
artifact in this region (Crenshaw, Bruegman, & Norman 1990). 
Weak C hi] emission at 3000 A may be present, but unfor- 
tunately it falls near a reseau and is thus partially obliterated. 
The ratio of the C iv/Lyman-a lines is nearly unity. The UV 
spectra are too noisy to detect emission-line variations at levels 
less than a factor of 3. 

To analyze the continuum variations, each individual spec- 
trum was corrected for effects of cosmic-ray hits and reseaux 
by linear interpolation between adjacent continuum regions 
and smoothed using a five-point boxcar. The SWP spectra 
were binned around 1600 A with a bin width of 600 A. The 
LWP binning was done at 2800 A with a bin width of 400 A. 

Figure 2 shows the UV long-term light curve at 1600 and 
2800 A. The 4400 A ( B band) light curve is also shown for 
comparison. The correlation between the frequency bands is 
apparent, although the UV points are sparsely sampled. The 
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TABLE 1 


Log of Observations 


MJD 

(1) 

Date 

(2) 

Observatory 

(3) 

Bandpass 

(4) 

8005 

1990 Apr 23 

UKIRT 

J, H, K, L 

8008 

1990 Apr 26 

Univ. Michigan 

4.8, 8 GHz 

8008 

1990 Apr 26 

Metsahovi 

14.5. 22. 37 GHz 

8310 

1991 Apr 22 

RHO, FBO 

B. V, R 

8315 

1991 Feb 27 

Metsahovi 

22, 37 GHz 

8315 

1991 Feb 27 

IUE 

1975-3300 A 

8316 

1991 Feb 28 

Univ. Michigan 

4.8 GHz 

8316 

1991 Feb 28 

IUE 

1175-2000 A 

8321 

1991 Mar 05 

Metsahovi 

14.5 GHz 

8321 

1991 Mar 05 

UKIRT 

J, H, K, L 

8322 

1991 Mar 06 

Univ. Michigan 

8 GHz 

8325 

1991 Mar 09 

Metsahovi 

22 GHz 

8326 

1991 Mar 10 

RHO. FBO 

B, V, R 

8326 

1991 Mar 10 

IUE 

1975-3300 A 

8339 

1991 Mar 23 

IUE 

1175-3300 A 

8342 

1991 Mar 26 

Metsahovi 

37 GHz 

8342 

1991 Mar 26 

RHO, FBO 

U, B, V, R 

8349 

1991 Apr 02 

Univ. Michigan 

4.8 GHz 

8351 

1991 Apr 04 

FBO 

R 

8353 

1991 Apr 06 

Metsahovi 

14.5 GHz 

8353 

1991 Apr 06 

Ginga 

1-20 keV 

8354 

1991 Apr 07 

Mestahovi 

22. 37 GHz 

8354 

1991 Apr 07 

IUE 

1175-3300 A 

8354 

1991 Apr 07 

IUE 

1175-3300 A 

8366 

1991 Apr 19 

UKIRT 

J, H, K, L 

8367 

1991 Apr 20 

Univ. Michigan 

4.8 GHz 

8367 

1991 Apr 20 

NRAO 

2 mm 

8367 

1991 Apr 20 

UKIRT 

0.8, 1.1, 1.3 mm 

8369 

1991 Apr 22 

Mestahovi 

37 GHz 

8369 

1991 Apr 22 

RHO, FBO 

U, B, V 

8369 

1991 Apr 22 

IUE 

1175-3300 A 

8370 

1991 Apr 23 

Metsahovi 

14.5, 22 GHz 

8371 

1992 Apr 24 

IUE 

1975-3300 A 

8393 

1991 May 16 

Univ. Michigan 

8 GHz 

8393 

1991 May 16 

Metsahovi 

22, 37 GHz 

8397 

1991 May 20 

Univ. Michigan 

4.8 GHz 

8397 

1991 May 20 

Metsahovi 

14.5 GHz 

8400 

1991 May 23 

UKIRT 

/ 

8400 

1991 May 23 

RHO, FBO 

B, V, R 

8402 

1991 May 25 

Ginga 

1-20 keV 

8404 

1991 May 27 

UKIRT 

/ 

8405 

1991 May 28 

Univ. Michigan 

4.8 GHz 

8405 

1991 May 28 

Metsahovi 

37 GHz 

8405 

1991 May 28 

RHO, FBO 

V, R 

8405 

1991 May 28 

IUE 

1975-3300 A 

8406 

1991 May 29 

Metsahovi 

14.5, 22 GHz 

8410 

1991 Jun 02 

Univ. Michigan 

8 GHz 

8414 

1991 Jun 06 

IUE 

1175-2000 A 

8417 

1991 Jun 09 

Metsahovi 

22, 37 GHz 

8417 

1991 Jun 09 

NRAO 

2 mm 

8417 

1991 Jun 09 

UKIRT 

0.8, 1.1, 1.3 mm 

8417 

1991 Jun 09 

RHO 

B 

8418 

1991 Jun 10 

UKIRT 

450 pm 

8424 

1991 Jun 16 

UKIRT 

J, H, K, L 

8447 

1991 Jul 09 

Metsahovi 

37 GHz 

8447 

1991 Jul 09 

UKIRT 

J , H, K, L 

8447 

1991 Jul 09 

UKIRT 

J, H, K, L 

8455 

1991 Jul 17 

NRAO 

2 mm 

8455 

1991 Jul 17 

UKIRT 

0.8, 1.1, 1.3 mm 


first two panels of Figure 3 show the 1UE UV observations 
made during the 1991 outburst. Figure 3 also shows the corre- 
spondence between the UV data and the B optical variations. 
The discrepant 1600 A (SWP) point which shows an increase in 
flux on 1991 March 23 (JD 2,448,338) is suspect because the 
signal-to-noise ratio in that observation was extremely low, 
and because a corresponding increase in brightness was not 
seen in the LWP camera. The IUE spectra are consistent with 


no variations in the UV line strengths throughout the decline 
of the outburst. 

2.3. Optical 

The optical data include 27 R-band CCD measurements 
from FBO and a large number of UBV photographic measure- 
ments from RHO. The RHO data were reduced as described in 
Webbet at. (1988). Additional CCD UBVRI observations were 
obtained with the 24 inch (0.6 m) telescope at Capilla Peak 
Observatory and the 1.9 m at Kitt Peak Observatory. The U-, 
B-, V-, R-, and /-band observations are plotted in individual 
panels of Figure 3. All magnitudes are dereddened following 
Seaton (1979), using a color excess of E(B— V) = 0.01 (Hewitt 
& Burbidge 1987). A comparison of the light curves shows that 
the outburst decline occurred concurrently in all seven fre- 
quency bands, although the entire UV-optical-IR region was 
well sampled for only four epochs during the flare decline. 

2.4. Infrared-Millimeter 

The IR data were obtained as part of an ongoing blazar 
monitoring program at the UKIRT. This program regularly 
observes a number of blazar sources with the UKIRT tele- 
scope at Mauna Kea (Robson 1992). The IR and 0.8-2.0 mm 
data were obtained as part of a long-term monitoring program 
of blazars using the UKIRT and the James Clerk Maxwell 
Telescope (JCMT). The infrared data were obtained using the 
UKIRT service program and used standard techniques of mea- 
surement and calibration. The JCMT data were obtained using 
the continuum photometer (UKT14) and were calibrated using 
planets and secondary calibration sources as described in 
Duncan et al. (1990). 

Observations made in the J, H, K, and L band (IR) are 
shown in Figure 4, while the flux variations at 800, 1100, and 
1300 p m and 2 mm are plotted in Figure 5. The J, H, K, and L 
bands showed a marked decrease in flux during the later 
phases of the optical outburst, followed by a small brightening. 
The 800, 1100, and 1300 pm fluxes showed a large increase 
corresponding to the outburst that was followed by a small dip 
in the flux curve. The variations are smoother at longer wave- 
lengths, although the sparse sampling prevented clear time- 
scales to be attributed to the flare rise. The 2 mm data also 
show a high flux level after the 1991 optical outburst, but no 
short timescale variations are seen in the flux curve. 

2.5. Radio 

Radio observations from the University of Michigan’s 26 m 
telescope were obtained at frequencies of 4.8, 8.0, and 14.5 
GHz. Additional radio observations at 22 and 37 GHz were 
obtained at the Metsahovi radio observatory in Finland 
(Terasranta et al. 1992). The radio flux curves are shown in 
Figure 6. It is worth pointing out that at frequencies 14.5, 22, 
and 37 GHz an upward light-curve inflection occurs nearly 
simultaneously with the onset of the optical flare. A similar 
upward inflection is seen at lower frequencies, 8 and 4.8 GHz, 
but at a later time, suggesting that the low-frequency emission 
emanates at larger distances from the central source than the 
higher frequency emission. An increase in flux at 8 GHz occurs 
approximately 200 days after the corresponding increase at 22 
GHz, and the increase at 4.8 GHz approximately 200 days 
after the 8 GHz flux begins to increase (Alter 1993). 

2.6. Multifrequency Spectra 

The data described in the previous sections have allowed us 
to construct 10 multifrequency spectra spanning at ~5 month 
baseline. The R-band light curve is plotted in Figure 7 with 
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Fig. 2. — Long-term UV and optical flux curves are plotted including all available archival 1UE data of 3C 345 and the RHO B-band observations. The B fluxes 
are plotted as dots, the LWP fluxes centered on 2800 A plotted as plus signs, and the SWP fluxes centered on 1600 A are represented by asterisks. This figure 
illustrates the correspondence between the optical flux variations and UV flux variations from about 1979 to the present outburst. 


plus signs along the bottom to indicate the approximate 
epochs at which these multifrequency spectra were compiled. 
The spectra were constructed from observations made over as 
short a time interval as possible (5-10 days), to maximize the 
frequency coverage while minimizing possible effects of short- 
timescale variability. Although 5 and 10 day bins are small 
compared with many previously published multifreqency 
spectra, significant variations do take place in some frequency 
bands within these time periods, thus introducing an addi- 
tional element of uncertainty into our analysis. 

We present the individual multifrequency spectra in Figure 
8. We have labeled the individual spectra “spl” through 
“ splO.” Table 2 is a compilation of the data that make up these 
spectra. The first row gives the spectrum identification (spl, 
etc.), while the second row indicates the date of the center 
around which all of the observations are clustered. The third 
row gives the half-width of the bin in days, and successive rows 
indicate the frequency bands in which data are collected and 
identified in column (1). The data are given as flux and 1 c 
errors and are in units of ergs cm' 2 s' 1 Hz' 1 . Two additional 
multifrequency spectra from Bregman et al. (1986) were 
analyzed along with the 10 spectra from this study. These addi- 
tional multifrequency spectra will be referred to as “sp81” and 
“ sp83.” A reference spectrum was created by adding these last 
two well-sampled spectra (sp81 and sp83); we consider this 
procedure reasonable because the brightness of 3C 345 was 
roughly constant during the 1981 and 1983 observations. The 
evolution of the multifrequency spectrum during the 1991 out- 
burst is illustrated in the succession of frames in Figure 8. The 
reference spectrum is plotted as small dots in each frame. The 


first spectrum, spl, was obtained prior to the outburst, while 
the other nine were made during the decline phase of the out- 
burst. 

Most of our data are consistent with a frequency- 
independent brightening of the source, the exception being the 
8 and 4.8 GHz fluxes. In particular, the onset of the increase in 
the flux at 14, 22, and 37 GHz occurs almost simultaneously 
with the onset of the optical flare. Although our millimeter and 
IR data do not cover this epoch, they are consistent with the 
onset of flaring having occurred during the same time interval. 
The fluxes at 8 and 4.8 GHz, while not varying synchronously 
with higher frequency emissions, do exhibit an increase which 
is readily interpreted and a delayed response to the primary 
flare. This behavior could be attributed to self-absorption at 
low frequencies, or could be indicative of origination at larger 
distances from the QSO central engine than the bulk of the 
emission. 

Examination of Figure 8 shows that the broad-band energy 
distribution derived from our observations has a slightly lower 
intensity and has a steeper IR-UV continuum than the refer- 
ence spectrum. The radio flux levels are similar to those of the 
reference spectrum. 

3. ANALYSIS AND MODELING 

The data collected during this outburst are ideal, for testing 
present theoretical models. Fitting homogeneous synchrotron- 
self-Compton (SSC) models (Urry 1984) to the data produced 
poor fits and unacceptably large Doppler factors (~ 100-150). 
In § 3.1 we examine the radio variability in detail, and in § 3.2 
the entire multifrequency spectral distribution is considered. 
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Fig. 3. — Individual UV (first two panels ) U, B. V, R, and I flux curves of 3C 345 during the decline of the 1991 outburst. The horizontal and vertical scales are the 
same in each frame, showing the relative activity in each band. 
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Fig. 3 — Continued 

3.1. Variability Modeling choice of indices T and p; however, the timescale for reaching 

The observed radio variations were compared with a simple ^e new e Quilibrium was different for the different electron 
model which follows the time evolution of the relativistic elec- energies owing to the energy-dependent character of synchro- 
tron distribution under conditions of continuous electron tron l° sses [dy/dt = by 2 and b = (4/3)a T c(B 2 /fto)]- When the 

injection, initiated at t = 0, in a volume of a given constant injected particles energies were folded through the synchrotron 

magnetic field. The resulting time-dependent synchrotron emissivity, a different characteristic relaxation timescale was 

emission was computed by folding the electron distribution obtained for the resulting synchrotron emission in different 

with the synchrotron emissivity in the <5-function approx- frequency bands. Under the conditions of continuous electron 

imation. At 1 = 0 the electron distribution function was injection and assuming the form of the electron spectrum at 

assumed to be a power law given by N e {y) — K e y~ p cm -3 t = 0 , the time-dependent electron distribution function is 

ergs - ', while starting at t = 0 a new population of electrons obtained by solving the electron kinetic equation. For bty < 1, 

was continuously injected with a power-law distribution of the distribution function is 

form Q e (y) = qy -r cm -3 s -1 ergs -1 , where T = p - 1 (r = 2 +J) 

in our simulations). The relationship between the indices of the N , ,\ — ir v -pm _ hrv) p ~ 2 4 - n — (1 — brv) r-1 "l 

injected and steady state electron distribution functions was * ’ ‘ h(T — 1) 

established so that after synchrotron losses the newly injected 

electrons developed a nonthermal power-law distribution with The distribution function is zero if bty > 1 by definition, 
spectral index identical to that of the electrons initially present When folded with the synchrotron emissivity, which can be 

in the given volume. As a result of these freshly injected elec- approximated by the <5- function e(v, v c ) cc <5(v - y 2 v c ), the 

trons, the steady electron distribution function at each energy above electron distribution function yields the synchrotron 

underwent an increase in amplitude. The change in amplitude emission at each frequency v, where v r = 4x 10 6 B Hz is the 

was the same for all electron energies because of the particular cyclotron frequency and B is in gauss. Because an electron of 



8000 8200 8400 8600 

JULIAN DAYS (-2,440,000) 

Fig. 4. — Same as Fig. 3, but for J, H, K, and L bands. The gap between 8000 and 8300 days prevents a determination of the rise time of the flare in those bands. 




576 


WEBB ET AL. 


Vol. 422 



Fig. 5. — Same as Fig. 3, but for 800, 1 100, 1300 /im and 2 mm bands. Again, the gap between 8140 and 8350 days allows only a lower limit estimate of a rise time. 


Lorentz factor y emits most of its energy at frequency y 2 v c and 
because the characteristic timescale for synchrotron losses are 
proportional to y -1 , one would expect in the context of this 
simple model that the variability timescale for each frequency 
would be inversely proportional to the square root of the fre- 
quency, i.e., t(v) oc v~ 1/j (see also Marscher & Gear 1985). The 
numerical modeling results fit this scaling law very well. 


The variability timescale for each radio-frequency band 
observed was estimated by fitting an exponential to the rising 
segment of the light curves. Unfortunately, our temporal 
coverage undersamples frequencies greater than 50 GHz (with 
the exception of optical frequencies). As a result, the variability 
timescales we obtain represent only upper limits to the real 
timescales. In the 2 mm and the 1300, 1 100, and 1300 ftm data, 



Fig. 6. — Same as Fig. 3, but for 37, 22, 14.5, 8.0, and 4.8 GHz data. The sampling in these bands allows an accurate determination of the rise times, and they are 
compared with theoretical models in § 3.1. 
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Fig. 7. — B-band light curve of 3C 345 during the 1991 outburst. The plus signs along the bottom of the plot indicate the approximate times when multifrequency 
spectra were compiled. The dashed line is intended merely to show the temporal order of the observations. The last two points are included to show a secondary 
outburst that occurred after our program has finished. 


there are large gaps which provide only upper limits on the 
variability timescales at those frequencies, but there is an indi- 
cation that these timescales tend to decrease toward higher 
frequencies also. We have calculated the relationship between 
the rise time and frequency from the 4.8, 8, 14, 22, and 37 GHz 
data and have plotted their logarithms in Figure 9. It is appar- 
ent that the relation obtained from these observations is in 
qualitative agreement with that expected within our simple 
model. Considering the simplifications in this model (uniform 
magnetic field, etc.), it is rather surprising that there is such 
good agreement with the data. Extrapolating the relationship 
thus obtained from the radio observations to the optical (using 
the scaling quoted before) yields a timescale of the order of 10 
days for the optical frequencies, a value consistent with some of 
the more rapid variations seen in the R-band optical data. 

3.2. Multifrequency Spectra Models 

In view of the success of the simple injection model, an effort 
was made to model the entire multifrequency spectrum in 
terms of an SSC model. Linear fits over restricted frequency 
ranges were made in order to determine the behavior of each 
part of the spectrum during the outburst. In the optical-IR 
region, the fits were subtracted from the IR-optical power-law 
component to obtain possible blue-bump components. If a 
discernible nonlinearity was present, we fitted a baseline to the 
data in that region instead of using a least-squares “best-fit” 
line to the data. The power laws obtained in the IR-optical-UV 
portion of the spectrum ranged from —1.27 to —2.7, with a 
mean of —1.58. Spectral indices derived in a similar fashion 
using the Bregman et al. (1986) multifrequency data give spec- 
tral slopes of —1.56 + 0.26 for the 1981 spectrum and 
— 1.49 + 0.07 for the 1983 spectrum. Brown et al. (1989) found 
the optical spectral index of 3C 345 to be between — 1.05 and 


- 1.74 with a mean of - 1.5, consistent with our observations. 
Although the statistics are poor, there seems to be no corre- 
lation of the spectral index with either flux level or time as the 
outburst progressed after subtraction of the blue-bump com- 
ponent. A comparison of the “ best-fit ” slopes (no bump com- 
ponent removed) yields essentially the same result. The 
complex multifrequency behavior exhibited by 3C 345 during 
the outburst decline suggests that more than one variable com- 
ponent is necessary in order to model the broad-band variabil- 
ity. 

VLBI studies of 3C 345 indicate the presence of a jet in 
3C 345 (Bartel et al. 1984; Tang, Ronnang, & Baath 1989). 
Bregman et al. (1986) and Sadun (1992) have found solid evi- 
dence for correlation between optical activity and super- 
luminal components in radio jets of some quasars. Due to the 
plausibility of relativistic motion in this source, we chose to 
model the broad-band continuum of 3C 345 in terms of a 
synchrotron-self-Compton jet (Konigl 1981 ; Hutter & Mufson 
1986) and the blue-bump component as a thermal accretion 
disk. 

The jet model assumes a power-law distribution of rela- 
tivistic electrons, N(e) = Kfr/r 0 )\ injected into a jet with 
opening half-angle <p, viewing angle 6, and no external pres- 
sure. The radius r is measured from the jet apex, and r 0 = 1.0 
pc. The electrons expand into a tangled magnetic field whose 
strength varies as a function of distance from the central source 
B = B,(r/r 0 ) m . Continuous reacceleration along the jet is 
assumed in order to offset radiative losses. 

One of the major problems with these types of jet models has 
been the large number of free parameters available to fit the 
data. Each individual model is poorly constrained by the 
observations, allowing the model to fit any spectrum. We have 
applied the model with the philosophy that once we obtain a fit 
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Fig. 8. — Multifrequency spectra throughout the 1991 outburst. The small dots represent the “reference spectrum" (see text), while the larger dots represent the observations at each epoch. Errors are 
smaller than symbol size unless explicitly shown. 










Fig. 8 — Continued 
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Data from Multi frequency Spectra 
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• The first line identifies the spectrum, the second line gives the date (see text), and the third line gives the bin half-width in days. 
b S = flux, dS = 1 terror, in units of ergs cm ' 2 s' 1 Hz' 1 . 
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Fig. 9. — The horizontal axis is log v, and the vertical axis is log a , where a 
parameterizes the frequency-dependent variability timescale r(v)acv - ’. The 
model predicts a = 0.5, which is depicted by the dashed line. The solid line 
represents a fit to the data with a = 0.54. 

to one of the well-sampled spectra, we then try to fit the other 
nine outburst spectra while adjusting the minimum number of 
free parameters, effectively reducing the number of uncon- 
strained parameters available to fit each individual spectrum. 

Conversion from flux to luminosity was done assuming a 
Friedmann cosmology with a q 0 of 0.0, a Hubble constant of 
100 km s' 1 Mpc“ \ and an emission-line redshift of 0.595. The 
continuum slope measured in each frequency band was used to 


transform each spectrum directly to the rest frame of 3C 345. 
The spectra sp5 and sp8 were chosen for the initial jet model 
fits because the Compton fluxes (X-ray observations) were 
available only for those spectra. The angles 9 and <f> are deter- 
mined by matching the optical and the X-ray flux levels in the 
models with those of our spectra. The turnover frequency v, 
and flux F t , the slope in the IR-optical region a sl , and the radio 
slopes a j2 were all constrained by fitting each spectrum. The 
Lorentz factor of the upper electron cutoff (T u ) was chosen so 
that the break flux and the optically thin flux either fit the 
entire IR-optical-UV region or fit the IR region while forming 
a baseline to the optical-UV flux. The parameters described 
above, along with the electron distribution index (a e ), are the 
input parameters for the jet models, and the numerical values 
of each spectrum are listed in rows 1-8 of Table 3. The second 
panel of Table 3 lists the calculated model parameters. The 
constants C x and C 2 are intermediate results, while K x and J9, 
are the normalizations for the electron distribution and mag- 
netic field as described previously. The values m and n are the 
exponents for the magnetic field and particle densities. R m is 
the inner radius of the emitting region in parsecs, and <5 is the 
Doppler boost factor. The third panel lists the mass and the 
accretion rate needed to account for the excess optical-UV flux 
in each spectrum which exhibited such an excess. The models 
calculated for each of the spectra are shown in Figure 10. The 
heavy dots are the data, the solid line represents the jet SSC 
spectrum, and the dotted line is the accretion disk flux. The 
turnover frequency was at lower frequencies, and the turnover 
flux was substantially less in spl than in the other spectra. 
Once the outburst began, the turnover frequency and flux 
increased and then remained effectively constant at log v 


TABLE 3 

Model Fit Parameters 


Spectrum 


Parameter sp81 sp83 spl sp2 sp3 sp4 sp5 sp6 sp7 sp8 sp9 splO 


A. Fitted SSC Jet Parameters’ 


6 0.11 0.06 0.14 0.14 0.14 0.14 0.14 0.14 0.14 0.14 0.14 0.14 

<t> 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 

V, 11.34 11.75 10.89 10.79 10.83 10.83 10.83 10.83 10.8 10.8 10.8 10.83 

F, -21.9 -22.1 -22.2 —21.8 -21.8 -21.8 -21.8 -21.8 -21.8 -21.8 -21.8 -21.8 

a. 0.59 0.69 0.59 0.59 0.59 0.59 0.59 0.59 0.59 0.59 0.59 0.59 

a„ 0.15 0.26 0.09 -0.39 -0.42 -0.42 -0.39 -0.39 -0.39 -0.39 -0.39 -0.39 

a, j 1.25 1.12 1.45 1.45 1.45 1.45 1.45 1.45 1.45 1.45 1.45 1.45 

r. 4.1 19 5.3 8.6 8.6 8.6 7.5 8.5 4.3 3.2 3.3 3.3 


B. Calculated SCC Jet Parameters b 


C, 1.2 51 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 

C 3 3.4 1.6 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 

m 0.87 0.96 0.85 1.12 1.15 1.15 1.12 1.12 1.12 1.12 1.12 1.12 

n 2.03 1.77 2.12 2.4 2.42 2.43 2.4 2.4 2.4 2.4 2.4 2.4 

JC, 1.33 4.47 1.1 1 1.4 1.37 1.2 1.2 1.2 1.2 1.2 1.2 

B t 0.83 0.56 1.16 2.6 2.4 2.5 2.45 2.45 2.6 2.6 2.6 2.6 

R m 1 0.34 2.72 7.3 6.6 7 6.8 6.8 7.4 6.4 7.4 7.4 

<5 6.76 16.52 4.53 4.53 4.53 4.53 4.53 4.53 4.53 4.53 4.53 4.53 


C. Fitted Accretion Disk Parameters' 


M 4.3 4.3 ... 4.5 4.3 4.5 4.3 4.3 ... 4.3 4.3 

M 3.3 2.2 ... 5.0 5.0 6.4 2.0 4.0 ... 1.0 4.1 


’ 6 and <5 are in radians; T is the Lorentz factor x 10’*. 

b C, is in units of 10 18 , C 2 is the units of 10~ u , JC, is in units of 10"*, R is in parsecs, and B is in gauss. 
* AJ is in units of 10~ * * * 6 * * 9 M 0 ,and Mis in M g yr~‘. 
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Fig. 10. — Model fits to the muftifrequency spectra, including the Bregman et al. (1986) spectra. The heavy dots are data points, the solid line indicates the Konigl (1981) SSC model, and the dotted line in 
the accretion disk component. 
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Fig. 10 — Continued 
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= 10.79 Hz and log F= —21.8 ergs cm -2 s" 1 Hz" 1 . The 
only parameter changed in generating jet models for sp2-spl0 
was the upper electron cutoff of the injected electron distribu- 
tion. The value of the upper electron cutoff energy changed 
from a Lorentz factor of 86,000 to one of 33,000 as the outburst 
subsided. Calculated values for the magnetic field strength, 
number density of electrons, and their radial dependencies all 
remaining fairly constant in the models as the outburst pro- 
gressed. 

3C 345 is among a subclass of blazar AGNs which exhibit 
the so-called big blue-bump component. In order to study the 
blue-bump emission, we subtracted the calculated best-fit jet 
model spectrum from the data to obtain residuals in the 
b)ue-UV spectral region. A statistically significant, positive 
residual is seen for six of seven cases where the blue-UV part of 
the spectrum is adequately sampled. We then fitted a thermal 
accretion disk spectrum to these residual spectra. The model 
involves thermal emission due to viscous dissipation in the 
steady state accretion disk which has been modified to include 
the effects of gravitational redshift and focusing (e.g., Sun & 
Malkan 1989); we have employed computational methods 
similar to those described in Czerny, Czerny, & Grindlay 
(1986). A Schwarzschild metric was used in the calculations, 
and an inclination of sin i = 0.5 was assumed. We fixed the 
inner disk radius at 3 R g and R 0ut /R in = 10 3 : 

Reasonable fits were obtained for several cases — sp2, sp3, 
sp83, and sp81 — with less satisfactory fits in other cases. There 
is no significant blue-bump excess in spl, sp7, and splO, so no 
fits were attempted. The fits to the multiepoch data were used 
to interpret the observed variability in the context of physical 
parameters. Initial fitting was used to estimate an “ average ” 
central black hole mass. We then fixed the mass to this value 
and repeated the fits, interpreting variations in the blue-bump 
luminosity as a varying accretion rate. 

Figure 11 depicts the residual blue-bump luminosity inte- 
grated over our fit to the data between log v = 14.5-16.5, 
versus the integrated underlying SSC jet luminosity. The upper 
two points (triangles) are from the Bregman et al. (1986) 
spectra. If there is any trend at all in our data, it is a slight 
increase in blue-bump luminosity with increasing SSC lumi- 
nosity. This suggests that if our interpretation of the blue bump 
as an accretion disk is correct, the accretion rate varies in the 
same sense as the SSC jet luminosity. Indeed, it is of interest 
that an excess is seen at all during outburst, since it is generally 
assumed that the synchrotron component is the dominant 
flaring component. We also note that the spectra for which no 
significant blue bump was seen extend over the full horizontal 
range of Figure 1 1. Table 3 lists the estimated parameter values 
and the calculated values in the context of the SSC jet/disk 
models. The spectra (including sp81 and sp83) are listed across 
the top, while the values of the various fitted and calculated 
parameters identified in column (1) are listed in vertical 
columns. Table 3A includes parameters constrained by direct 
model fitting to the data. The calculated physical parameters 
such as the magnetic field and minimum radius are given in 
Table 3B. Table 3C gives the mass and accretion rate derived 
when possible. It should be noted that the 1981 and 1993 fitted 
and calculated parameters were all very different from those 
derived for the 1991 flare. 


3C 345 - Blue Bump Residuol vs SSC Jet luminosity 



F 10 . 11. — “Blue-bump" luminosity vs. Ihe SSC-je! luminosity. A mild 
correlation exists between the luminosities for the 1991 outburst data. The two 
upper points are the sp81 and sp83 spectra. 


4. CONCLUSIONS 

We have undertaken a program to study the multifrequency 
variability of a blazar during an energetic outburst. From our 
observations of the 1991 3C 345 outburst, we have constructed 
10 nearly simultaneous multifrequency spectra, some of which 
cover a range of nearly 10 decades in frequency. Comparison of 
the relationship between frequency and rise time in a simple 
electron injection model and the rise times seen in the radio 
flux curves show that the radio data is consistent with the 
simple injection model. The multifrequency spectral evolution 
was then modeled in terms of a SSC jet in which the upper 
cutoff of the relativistic electron distribution changes with time, 
and an accretion disk in which the accretion rate also varies 
with time. Modeling several multifrequency spectra of the same 
source during an outburst allows us effectively to decrease the 
number of free parameters, which is a common problem with 
complicated jet models. Comparison of the luminosity in the 
jet component with that of the blue-bump component shows a 
marginal linear correlation, suggesting that the accretion rate 
increases as the jet luminosity increases. 

Further separation of the jet and disk components would be 
expedited by obtaining simultaneous optical polarization mea- 
surements during an outburst. Data of this type, with sufficient 
frequency and temporal sampling, could provide stringent tests 
of more physical jet models, i.e., those incorporating shocks 
and using numerical simulations instead of analytical approx- 
imations (such as reacceleration). 
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